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ABSTRACT 


Within the lower non-dissipative atmosphere, the eigenfunctions of 
the atmosphere are the Hough-functions which describe the tidal wave motions. 
Under the influence of dissipation mechanisms like heat conduction, viscosity 
and neutral-ion drag these Hough -functions become modified at thermospheric 
heights and are shown to degenerate ultimately into spherical harmonics at 
heights above 200 km. The implications of this degeneration on the latitudinal 
structure of the tidal modes and on the vertical wave propagation are dis- 
cussed. It is suggested in particular that the semi-diurnal mode becomes 
quasitrapped at altitudes above 160 km thus explaining the apparent absence 
of semi-diurnal variations in the dynamics of the upper thermosphere. 


A common way to treat thermosphere dynamics is to separate the varia- 
bles and this is only possible under the assumption that perturbation theory 
is a good approximation. One asks for solutions of the form 

F(z) 0 (6) A (A,t) (1) 

where F is a function of height z, 0 is a function of co -latitude <5 and A 
is a function of longitude A and universal time t. In the case of regular 
tidal and planetary waves, A is a simple exponential function 

A = exp [jmA + jwt] (2) 

where m is an integer and o) is the angular frequency of the wave. It remains 
then to determine the two functions F and 0 In (1) . 

Since the earth’s atmosphere is a cavity, it behaves like a spherical 
wave guide in which only discrete geometrical structures of the waves can 
exist. The eigenfunctions that describe these wave modes are the Hough-func- 
tions within the lower atmosphere [1] and are identical with the function 0 
in (1) in this height region. Hough- functions can be represented by a series 
of spherical harmonics and generally have a rather complicated latitudinal 
structure due to the strong influence of the Coriolis force on the wave 
propagation. However, within the thermosphere, dissipation effects like heat 
conduction, viscosity and neutral -ion drag become increasingly important with 
increasing heights and reduce the influence of the Coriolis force. The 
Hough- functions therefore change their latitudinal structures with altitude 
between about 100 and 200 km and ultimately degenerate into the spherical 
harmonics above about 200 km [2] : 


above 200 km 




(3) 


The indices n and m of tne function 0 are identical with the wave domain 
numbers n and m of the corresponding spherical function p 1 ^ , and 

f * w/2ft (4) 

is the Coriolis parameter (ft is the angular frequency of the earth's rotation). 

Figure 1 shows as an example the pressure and the horizontal wind compnn- 

i > 1/2 

ents of the fundamental symmetric diurnal tidal wave mode of gravity 

wave type. In Figure la the pressure function of this mode is plotted versus 
co-latitude for four different heights between 100 and 300 km. At 100 km 
altitude a non-dissipative atmosphere has been assumed. Here the Hough-function 
is the trapped (l,-l)-m'ode of Kato [3] or the (1,-2) -mode of Lindzen [4]. The 
dash -dotted line of the pressure function at 100 km gives the exact structure 
of this mode while the solid line represents an approximation introduced to 
simplify the analytic solutions. The dashed lines give the phase of the (gen- 
erally complex) pressure function which is zero at 100 km indicating that 0 is 
a real function there. 

Between 100 and 200 km, the structure of the pressure function changes 

until finally it is transformed into the spherical harmonic which describes 

there the thermospheric diurnal pressure bulge of the Jacchia -model [5] . 

In Figure lb and lc, the meridional and zonal wind components (blowing 

1 1/2 

from the north and from the west, respectively) of the diurnal mode 0^* 

are plotted versus co- latitude in the same manner as in the case of the pressure 

function in Figure la. Here the zonal wind at 100 km altitude and below changes 


its direction at about 30° latitude which is indicated by a phase jump o£ 

12 hours. This wind field is mainly responsible for the generation of the 
geomagnetic Sq current at E- layer heights [3], [6], [7] . 

During its transition within the height range between 100 and 200 km, 
this reversal of the zonal wind direction disappears, and the horizontal wind 
field above 200 km is essentially identical with the wind field derived by 
Geisler [8] and by Kohl and King [9] from the Jacchia-model [5]. 

The transition region between 100 and 200 km of this wave mode and likewise 
of all the other wave modes , too , imposes a problem in connection with our 
basic assumption that the variables in Eq. (1) can be separated. As we note 
from Figure 1, within this height range the latitude functions 0 depend not 
only on co- latitude but also on height which contradicts this assumption. 
Apparently, a coupling between the various wave modes is to be expected in that 
height range. The degree of coupling can only be determined by an exact solution 
of the system of partial differential equations without separation of the varia- 
bles. If we nevertheless do such a separation as in Figure 1 and consider the 
different wave modes as if they were decoupled from each other we in fact make 
a kind of ray optics approximation well known from the electromagnetic wave 
propagation within the ionosphere. 

A coupling between the various wave modes takes place although perturbation 
approximation has been assumed. Taking into account the non-linear terms in the 
hydrodynamic equations, gives rise to an additional coupling between the modes 
and complicates the problem furthermore [10] , [11] . 

Assuming the separation of the variables as a sufficient approximation, 
one obtains separation constants for the various wave modes called the equivalent 
depths. These equivalent depths are a measure for the eigenvalues of the height 


3 


function F in Eg. (1). Taking into account heat conduction, viscosity and 
neutral -ion drag, gives 8 complex eigenvalues corresponding to 8 character- 
istic solutions or characteristic waves of the system of the 8 complex coupled 
ordinary differential equations of first order for the height function F of 
each wave mode. Within a thin strip at a certain height in which the physical 
parameters of the thermosphere can be considered as constant, the character- 
istic waves are decoupled from each other, and each characteristic wave has 
the height structure 

exp[-(jou +3^-13 z/ 2H q ] (5) 


where and 3^ are real and imaginary part, respectively, of the eigenvalues 


• th 


of the i characteristic wave of the wave mode (n,m,f) , and H is the pressure 


scale height within this strip. 

Within the lower atmosphere, only characteristic waves of gravity wave 
type can exist. They are either propagation modes (a 0; 3^ * 0) or trapped 
modes (c^ = 0; 3^ ^ 0) . The (1,-1) -mode is such a trapped mode with [3^ | > X 
Within the thermosphere, heat conduction waves and two kinds of viscosity waves 
can exist besides the gravity waves. Eacu wave type consists of an upgoing 
wave (3- >0) and a downgoing wave (3* < 0). Above about 200 km, all character 


istic waves are quasi -trapped 0; |3^ | > 1) . 


In Figure 2, we plotted versus altitude the attenuation factors (3 -1 and 
the propagation factors a of the upgoing gravity wave type of three tidal 
wave modes - the fundamental symmetric diurnal mode Sl* 1 ^ 2 or the (1,-1) -mode 


of Kato [3], the fundamental symmetric semi-diurnal mode 0 2,1 or the (2,2) 


-mode according to Chapman and Lindzen [1] and the fundamental antisymmetric 
diurnal mode 0 | ,1 ^ 2 * The ©j* 1 ' 2 -mode being a trapped mode within the lower 


. Wii , , a fi t ' " ■ _ ■ 
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atmosphere remains a quasi-trapped mode at thermospheric heights with de- 
creasing attenuation factor above 140 km. The semidiurnal ©I* 1 ” m °de ^* s a 
propagation mode at 100 km (0 ~ 0) . However, with increasing height, the 
attenuation factor 3 decreases rapidly (the dashed line in Figure 2a which 
indicates negative values of 3 - 1) and becomes greater than one at 160 km. 
Above this height, the G^’ 1 -mode behaves like a quasi-trapped mode. The 
antisymmetric G^ 1 ^ 2 -mode which describes the annual modulation of the diurnal 
variation of pressure and winds has the eigenvalue (3=1; a = 0) witlun the 
lower atmosphere. This mode too changes its character into a quasi-trapped 
mode above 100 km. 

The propagation factors a which are a measure for the vertical phase 
velocities increase with altitude between 100 and 150 km and then decrease 
reaching assymptotic values above about 300 km. At each height one can define 
a vertical wave length according to 

X z = 4 it H q / g, * (6) 

The upper abscissa in Figure 2b is scaled according to X /H and shows 
that at 150 km the tidal waves have minimum vertical wave lengths. 

Figure 2 shows only the one most important upgoing characteristic gravity 
wave for the three respective tidal modes. This characteristic wave together 
with the seven others can be excited in any height by a heat source component 
corresponding to the respective tidal mode. The eight characteristic waves 
propagate up and down and become coupled with each other at any height within 
a realistic thermosphere. 

If one considers a heat input limited within a certain height region 
between z q and z^, the boundary conditions are the radiation conditions which 


.imply zero wave energy flux from outside through the boundaries (z , 2 ^) of 


the model. These conditions are fulfilled by setting the incoming character' 
istic waves from above (that are waves with 8 .< 0) equal to zero at the 

JU 

upper boundary z^ and by setting waves from below (8^ > 0) equal to zero at 
the lower boundary a . 

Such procedure can be handled without great mathematical difficulty 


in numerical calculations if one neglects the viscosity waves [1,2] , Other- 
wise, one must use the conditions of zero flux of mass, momentum and energy 
through the upper boundary as an alternative upper boundary condition [13] , 
[14] , [15] by which one looses however the advantage of the discriminating 
various characteristic waves. 
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